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^ Abstract 

' We reanalyse the potential of the LHC to discover new Z' gauge bosons and to 

^ ■ discriminate between various theoretical models. Using a fast LHC detector simu- 

lation, we have investigated how well the characteristics of Z' bosons from different 
models can be measured. For this analysis we have combined the information com- 
ing from the cross section measurement, which provides also the Z' mass and total 
width, the forward-backward charge asymmetries on- and off-peak, and the Z' ra- 
pidity distribution, which is sensitive to its uu and dd couplings. We confirm that 
new Z' bosons can be observed in the process pp — > Z' — > £ + £~, up to masses of 
about 5 TeV for an integrated luminosity of 100 fb _1 . The off- and on-resonance 
peak forward- backward charge asymmetries show that interesting statistical 
accuracies can be obtained up to Z' masses of the order of 2 TeV. We then show 
how the different experimental observables allow for a diagnosis of the Z' boson and 
the distinction between the various considered models. 
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1. Introduction 



Although the Standard Model (SM) of the electroweak and strong interactions describes 
nearly all experimental data available today [1], it is widely believed that it is not the 
ultimate theory. Grand Unified Theories (GUTs), eventually supplemented by Supersym- 
metry to achieve a successful unification of the three gauge coupling constants at the high 
scale, are prime candidates for the physics beyond the SM. Many of these GUTs, including 
superstring and left-right-symmetric models, predict the existence of new neutral gauge 
bosons, which might be light enough to be accessible at current and/or future colliders 1 ; 
for reviews see Ref. [3]. New vector bosons also appear in models of dynamical symmetry 
breaking [4] and recently, "little Higgs" models have been proposed to solve the hierarchy 
problem of the SM [5]: they have large gauge group structures and therefore predict a 
plethora of new gauge bosons with masses in the TeV range. 

The search for these Z' particles is an important aspect of the experimental physics 
program of future high-energy colliders. Present limits from direct production at the 
Tevatron and virtual effects at LEP, through interference or mixing with the Z boson, 
imply that new Z' bosons are rather heavy and mix very little with the Z boson. Depend- 
ing on the considered theoretical models, Z 1 masses of the order of 500 to 800 GeV and 
Z-Z 1 mixing angles at the level of a few per- mile are excluded 2 [7] . A Z 1 boson, if lighter 
than about 1 TeV, could be discovered at Run II of the Tevatron [8] in the Drell-Yan 
process pp — > Z' — > £ + £~, with £ — e,/x [9]. Detailed theoretical [8] and experimental 
[10, 11, 12] analyses have shown that the discovery potential of the LHC experiments is 
about 5 TeV, using the process pp — > Z' — > £ + £~ . Future e + e~ colliders with high cm. 
energies and longitudinally polarized beams could indicate the existence of Z 1 bosons via 
its interference effects, with masses up to about 6 x y/s [8, 13]. 

After the discovery of a Z 1 boson, some diagnosis of its coupling needs to be done in 
order to identify the correct theoretical frame. For this purpose, and since a long time, 
the forward-backward charge asymmetry for leptons Ap B has been advocated as being a 
powerful tool [14]; the most direct method to actually measure Ap B at the LHC has been 
described in [15]. In addition to the information from the total Z' cross section, it has 
been argued that the measurement of ratios of Z ! cross sections in different rapidity bins 
might provide some information about the Z 1 couplings to up and down quarks [16]. 

Following the arguments given in [17], we advocate that the Z 1 cross section should be 
measured relative to the number of produced Z bosons for the same lepton final states. 
Using this approach, many systematic uncertainties due to theoretical and experimental 
uncertainties will cancel, and the relative Z' jZ cross section ratio might be measured 

1 For example the breaking at the supersymmetry-breaking scale, i.e. at a scale around the TeV, of 
an extra U(l)/ group to which a Z' boson is associated, might solve the so-called problem, which 
notoriously appears in the Minimal Supcrsymmctric extension of the SM [2] . 

2 In contrast, some experimental data on atomic parity violation and deep inelastic neutrino-nucleon 
scattering, although controversial and of small statistical significance [see Ref. [1] for instance], can be 
explained by the presence of a Z' boson [6] . 
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and calculated with an accuracy of about 1%. Furthermore, the method should also 
lead to precise relative parton distribution functions for u and d quarks, as well as for the 
corresponding sea quarks and antiquarks. Thus, we can go beyond the previously proposed 
procedure to analyse the Z' rapidity distribution [16], by performing a fit. The fit uses the 
predicted rapidity spectra as calculated with uu and dd, as well as the contribution of the 
sea, for the mass region of interest, which is directly related to xi, X2 of the corresponding 
quarks and antiquarks in the proton. 

While numerous theoretical and experimentally motivated Z' studies have already 
been performed, the combination of all sensitive LHC variables, as described above, has 
not been done so far; the work described in this paper will thus fill a gap. We will 
perform the studies using the PYTHIA program [18] and a fast LHC detector simulation. 
We first update previous studies using, the latest parton distribution functions [21], and 
extend them in two directions. First, following the method proposed in [15], the forward- 
backward charge asymmetries, on and off the Z 1 resonance peak, are analysed together 
with the cross section in order to differentiate between the different models 3 . Second, 
we show that a direct fit of the rapidity distribution allows for additional information 
and would be useful to disentangle between Z' bosons from various models through their 
different couplings to up-type and down-type quarks. 

The rest of the discussion will be organized as follows. In the next section, we define the 
theoretical framework in which our analysis will be performed. In section 3, we describe 
the relevant observables that can be measured at the LHC, namely the dilepton cross 
section times the Z 1 total width, the on-peak and off-peak forward-backward asymmetries 
and the rapidity distribution, and the simulation tools which we will use in our study. In 
section 4, we analyse the resolving power of these observables. 

2. The considered Z' models 

To simplify the discussion, we will focus in this paper on two effective theories of well 
motivated models that lead to an extra gauge boson: 

1) An effective SU(2)l x U(1)y x U(1)y' model, which originates from the breaking of 
the exceptional group E 6 , which is general enough to include many interesting possibilities. 
Indeed, in the breaking of this group down to the SM symmetry, two additional neutral 
gauge bosons could appear. For simplicity we assume that only the lightest Z' can be 
produced at the LHC. It is defined as 

Z' = ^cos/3 + Z;sin/3 (1) 

and can be parametrized in terms of the hypercharges of the two groups U(l)^ and U(l) x 
which are involved in the breaking chain: 

E 6 -> SO(10) x U(l)^ -> SU(5) x U(l) x x U(l)^ -> SU(3) C x SU(2) L x U(1) Y x U(l)y 

3 Recently, the off-peak forward-backward asymmetry has also been used in Ref. [19] to study Kaluza- 
Klein excitations of gauge bosons. 
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The values (3 = and (3 = tt/2 would correspond, respectively, to pure Z' x and Z'^ 
bosons, while the value (3 = arctan(— \Jb/3) would correspond to a Z' v boson originating 
from the direct breaking of E 6 to a rank-5 group in superstrings inspired models. 

2) Left-right (LR) models, based on the symmetry group SU(2) R x SU(2) L x U(1)b-l, 
where B and L are the baryon and lepton numbers. Even though we investigate only the Z' 
in this paper, it should be recalled that new charged vector bosons, potentially observable 
at the LHC, also appear in these models. The most general neutral boson Z' LR will couple 
to a linear combination of the right-handed and B—L currents: 



Jlr = (XlrJsr - (l/2a LR )J£_ L with a LR = yf (c 2 w g 2 R / ' s 2 w g 2 L ) - 1 , 



(2) 
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where gL—^/sw an d 9r are the SU(2)l and SU(2)r coupling constants with s\ 
sin 2 9 W . The parameter a LR is restricted to lie in the range ^2/3 < a^ R < \/2: the upper 
bound corresponds to a LR-symmetric model with g R = gi, abbreviated in the following 
as LR, while the lower bound corresponds to the x model discussed in scenario (1), since 
SO(10) can lead to both SU(5) x U(l) and SU(2) R x SU(2) L x U(l) breaking patterns. 

In order to achieve a complete comparison, we will also discuss the non-realistic case 
of a sequential boson Z$ M , which has the same fermion couplings as the SM Z boson, 
as well as a Z' boson, denoted by Z' d , with vanishing axial and vectorial couplings to u 
quarks and which, in E Q models, corresponds to the choice cos/9 = v/o/8. 



The left- and right-handed couplings of the Z' boson to fermions, defined as: 
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are given in Table 1 for the first-generation fermions in the two scenarios. The mixing 
between the Z and Z' bosons is very small [7] and will be neglected in our discussion. 
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Table 1: Left- and right-handed couplings of the Z' boson to the SM fermions with the 
notation of the first generation in the Eq (left panels) and LR (right panels) models. 

The Z' partial decay width into a massless fermion-antifermion pair reads 

aM z > 



(iff + (9 f / f 



(A) 



with N c the colour factor and the electromagnetic coupling constant to be evaluated at 
the scale Mz> leading to a ~ 1/128. In the absence of any exotic decay channel, the 
branching fractions for decays into the first-generation leptons and quarks are shown in 
Fig. 1 for Eq and LR models as functions of cos/3 and am, respectively. As can be seen, 
the decay fractions into £ + £~ pairs are rather small, varying between 6.6% and 3.4% for Eq 
models and 6.6% and 2.3% for LR models; in the latter case the decay branching fraction 
is largest for the symmetric case = g^ and smallest for a^R — \Pl. The Z' total decay 
width, normalized to Mz', is also shown in Fig. 1: it is largest when cos/3 = ±1 in E e 
models and ~ \[2 in LR ones. The Z' bosons that we will consider here are thus 
narrow resonances, since their total decay width does not exceed 2% of their masses 4 . 
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Figure 1: The branching ratios of the decays Z' — > ff in E Q models as a function of 
cos/3 (left) and in LR models as a function of oilr (right). The total Z 1 decay widths, 
normalized to 10 /Mz>, are also shown. 

In the limit of negligible fermion masses, the differential cross section for the subprocess 
qq — > £ + £~ , with respect to 9* defined as the angle between the initial quark q and the final 
lepton £~ in the Z' rest frame, is given by [s = Mf e is the cm. energy of the subprocess] 



-{qq^j,Z, Z' 



1 wa z 



dcos^* v " " ' ' 9 2s 

where the charges Qi and Q3 are given by [13] 



;i + cos 2 0*)Q 1 + 2cos6*Qi 



Qi/3= \Qll\ 2 + \Qrr\ 2 ±\Qrl\ 2 ±\Qlr\ 2 /4 



(5) 



(6) 



4 Note however that non-standard decays, such as decays into supersymmetric particles and/or decays 
into exotic fermions, are possible; if kinematically allowed, they can increase the total decay width and 
hence decrease the Z' — > £ + £~ branching ratios. In the case of the Eq model for instance, the fermions 
belong to a representation of dimension 27 which contains 12 new heavy states per generation, and if 
they are light enough, the total decay width of the Z' is then simply Y z> — 2.5Mz*% independently of 
the angle (3 [14]. These exotic fermions, however, should also be observed at future colliders; see e.g. [20]. 
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In terms of the left- and right-handed couplings of the Z' boson defined previously, and of 
those of the Z boson = I^ — Q^s^, = —QfSyy] and the photon [g{ 7 = — Qf] 
with Qf the electric charge and l[ L the left-handed weak isospin, the helicity amplitudes 
Qij with i,j — L,R for a given initial qq state read 

a qZ (/ z S n qZ ' n lz ' S 

qi. = ri h + 9 -^l± - + 9i 9 -i - (7) 

WtJ 9t 9 3 s* w c 2 w S-MI + iT z M z c* w S - Ml + iT z ,M z , ' 1 ] 

To obtain the total hadronic cross section 5 and forward-backward asymmetries, we must 
sum over the contributing quarks and fold with the parton luminosities. 

A few points are worth recalling concerning the forward-backward asymmetry in E$ 
models [14]: (i) since the up-type quarks have no axial couplings to the Z' boson, Q\ = 0, 
they do not contribute to Ap B on the Z' peak; (ii) the asymmetry completely vanishes 
for three (5 values: (5 = arctan(— ^3/5) and (3 = ±7r/2, where the left- and right-handed 
Z' couplings of both rf-quarks and charged leptons are equal; (in) off the Z' resonance, 
there is always an asymmetry that is generated by the Z boson couplings. 



3. Observables sensitive to Z' properties 

The LHC discovery potential for a Z' as a mass peak above a small background in the 
reaction pp — > Z' — > with £ — e, fx, is well known. The required luminosity to 

discover a Z' basically depends only on its cross section, and therefore on its mass and 
couplings. Experimental effects due to mass resolution, assuming the design parameters of 
ATLAS or CMS [11, 12], are known to result in an only minor reduction of the sensitivity. 

Once a Z 1 boson is observed at the LHC, we will obviously measure its mass, its total 
width and cross section. Furthermore, forward-backward charge asymmetries on and 
off the Z' resonance provide additional information about its couplings and interference 
effects with the Z boson and the photon. In addition one can include the analysis of the 
Z' rapidity distribution, which is sensitive to the Z' couplings to uu and dd quarks. Such 
future measurements can be performed as follows at the LHC: 

The total decay width of the Z' is obtained from a fit to the invariant mass distribu- 
tion of the reconstructed dilepton system using a non-relativistic Breit-Wigner function: 
a /[(M! - Ml,) 2 + <ii] with a x = T|,M|,. 

The Z' cross section times leptonic branching ratio is calculated from the number 
of reconstructed dilepton events lying within ±3r around the observed peak. The 3r 
interval used to define the cross section is arbitrary; however, if varied from 2r to 5r, the 
cross section increases only between 5 and 10% for different Z' models and masses 6 . 

5 A if- factor of the order of JsTdy ~ 1.4 [22] for the production cross section can be also included. 

6 As noted previously, both the total width and the cross section times the leptonic branching ratio 
can be altered if exotic decays of the Z' boson are present. However, this dependence disappears in the 
product, and it is this quantity that should be used in discriminating models independently of the decays. 
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The leptonic forward-backward charge asymmetry A l FB is defined from the lepton an- 
gular distribution with respect to the quark direction in the centre-of-mass frame, as: 



da 



oc 



1 + cos 2 r) + Aiocosfl* . 



d cos 9* 8 



FB 



(8) 



As discussed in Ref. [15], the lepton angle 9* in the dilepton cm. frame can be calculated 
using the measured four momenta of the dilepton system. A l FB can then be determined 
with an unbinned maximum likelihood fit to the cos 9* distribution. Figure 2a shows the 
asymmetry, assuming that the quark direction is known, as a function of the dilepton mass 
for a Z' and a Z^ M boson, assuming a mass Mz< = 1.5 TeV. A FB varies strongly with 
the dilepton mass and is very different in the two models. Unfortunately, A FB cannot be 
measured directly in a proton-proton collider, as the original quark direction is not known. 
However, it can be extracted from the kinematics of the dilepton system, as was shown in 
detail in [15]. The method is based on the different x spectra of the quarks and antiquarks 
in the proton, which allows to approximate the quark direction with the boost direction 
of the ii system with respect to the beam axis (the z axis). Consequently, the probability 
to assign the correct quark direction increases for larger rapidities of the dilepton system 
and somewhat cleaner and more significant measurements can be performed, as shown in 
Fig. 2b. A purer, though smaller, signal sample can thus be obtained by introducing a 
rapidity cut. For the following studies we will require \Yu\ > 0.8. 
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Fi gure 2: Ap B as a function of the dilepton mass: the predicted asymmetry for a Z' and 
a Z' SM together with the SM prediction (a) and in different dilepton intervals (b). 



The Z' rapidity distribution allows us to obtain the fraction of Z' bosons produced 
from uu and dd initial states. Assuming that the W ± and Z boson rapidity distributions 
have been measured in detail, as discussed in [17], relative parton distribution functions 
for u and d quarks, as well as for the corresponding sea quarks and antiquarks are well 
known. Thus, the rapidity spectra can be calculated separately for uu and dd, as well 
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as for sea quark antiquark annihilation, and for the mass region of interest. Using these 
distributions, a fit can be performed to the Z 1 rapidity distribution, which allows to 
obtain the corresponding fractions of the Z 1 boson produced from uu, dd as well as for 
sea quark-antiquark annihilation 7 . 

In the present analysis, PYTHIA events of the type pp — > 7, Z, Z' — > ee, fifi were simu- 
lated at a centre-of-mass energy of 14 TeV, and for the Z' models discussed in section 2. 
The Z' masses were varied from 1 TeV up to 5 TeV. These events were analysed, using 
simple acceptance cuts following the design criteria of ATLAS and CMS. Following the 
results from previous studies and the expected excellent detector resolutions, the obtained 
values are known to be rather insensitive to measurement errors, especially for the e + e~ 
final states. We therefore do not include any resolution for the current study. In detail, 
the following basic event selection criteria were used: 

• The transverse momenta of the leptons, p l T , should be at least 20 GeV. 

• The pseudorapidity \rj\ of each lepton should be smaller than 2.5. 

• The leptons should be isolated, requiring that the lepton carries at least 95% of the 
total transverse energy found in a cone of size of 0.5 around the lepton. 

• There should be exactly two isolated leptons with opposite charge in each event. 

• The two leptons should be back to back in the plane transverse to the beam direction, 
so that the opening angle between them was larger than 160°. 

Figure 3 shows the expected number of events for a luminosity of 100 fb -1 . The SM 
background relative to the signal cross section is found to be essentially negligible for the 
considered Z' models. We thus reconfirm the known Z' boson LHC discovery potential, 
to reach masses up to about 5 TeV for a luminosity of 100 fb -1 [8]. 

Figure 4a shows the invariant mass distribution for the dilepton system, as expected 
for different models with M z > fixed to 1.5 TeV and for the SM using a luminosity of 100 
fb -1 . For all Z' models, huge peaks, corresponding to 3000-6000 signal events, are found 
above a small background. The cross sections for Z 1 bosons in the various models are also 
strongly varying. 

In addition, very distinct observable forward-backward charge asymmetries are ex- 
pected as a function of the dilepton mass and for the different Z 1 models, as shown in 
Fig. 4b. In order to get an impression of how an experimental signal with statistical 
fluctuations would look like, the measurable forward-backward asymmetry in the Z' v case 
has been generated with the number of events corresponding to 100 fb -1 , as shown in 
Fig. 4b. We find that additional and complementary informations is also obtained from 
Ap B measured in the interference region. To quantify the study for a Z' mass of 1.5 TeV, 

7 Following this procedure, and having very large statistics at hand, it would be imaginable even to 
measure also the forward-backward charge asymmetries separately for u and d quarks. Charge asym- 
metries for different Z' rapidity intervals would have to be measured and, with the knowledge of the 
corresponding uu and dd fractions from the entire rapidity distribution, the corresponding u and d asym- 
metries could eventually be disentangled. However, a quick analysis of the potential sensitivity indicates 
that an interesting statistical sensitivity would require a luminosity of at least 1000 fb -1 . 
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6 Accepted events for 100 fb" 1 
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Figure 3: The LHC discovery reach: the number of events fulfilling the selection cuts for 
a luminosity of 100 fir 1 , for the different Z' models and in the SM case. 

"on-peak events" are counted if the dilepton mass is found in the interval 1.45 TeV 

< Mm < 1.55 TeV. The "interference region" is defined accordingly and satisfy 1 TeV 

< M u < 1-45 TeV. 

Finally, the rapidity distribution is analysed. Figure 5a shows the normalized distri- 
butions for a Z' with a mass of 1.5 TeV produced from uu, dd and sea-antisea quark 
annihilation. Especially the Z' rapidity distribution from uu annihilation appears to be 
significantly different from the other two distributions. Figure 5b shows the expected ra- 
pidity distribution for the Z' v model. A particular Z' rapidity distribution is fitted using 
a linear combination of the three pure quark-antiquark rapidity distributions shown in 
Fig. 5b. The fit output gives the uu, dd and sea quarks fraction in the sample. This will 
thus reveal how the Z' couples to different quark flavours in a particular model. 

In order to demonstrate the analysis power of this method, we also show the rapidity 
distribution in the case of the Z^ boson, which has equal couplings to up-type and down- 
type quarks. As can be qualitatively expected from the distributions shown in Fig. 5a, 
the used fitting procedure provides very accurate results for the known generated fraction 
R U u of uw/all, while some correlations between dd and the sea-antisea Z' production, 
which limits the accuracy of the measurement for the dd fractions. For example, for the 
Z' v model, the generated event fractions from uu, dd and sea-antisea quarks are 0.71, 
0.26 and 0.03 respectively. The corresponding numbers from the fit and 100 fb -1 are 
0.71±0.07, 0.29±0.08 and 0.01±0.02. 

Table 2 shows the value of the cross section times the total decay width, the forward- 
backward charge asymmetry for the on-peak and interference regions as defined above, 
and the ratio of Z' events produced from uu annihilation as obtained from the fit to the 
Z' rapidity distribution. 
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Dilepton invariant mass spectrum Forward backward asymmetry measurement 
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Figure 4: The dilepton invariant mass spectrum (a) and A e FB (b) as a function of Mu for 
four Z' models. For the forward-backward charge asymmetry, the rapidity of the dilepton 
system is required to be larger than 0. 8. A simulation of the statistical errors, including 
random fluctuations of the Z' v model and with errors corresponding to a luminosity of 
100 fb^ 1 has been included in (b). 




Figure 5: The normalized rapidity distribution of Z' with a mass o/1.5±0.05 TeV produced 
from the different types of quarks (a). The observable rapidity distribution for two different 
Z' models is shown in (b ), including the fit results that determine the different types of qq 
fractions. 
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Table 2: The values of the four basic observables, the signal cross section, multiplied by 
the total width, the forward-backward charge asymmetry on- and off-peak, and the ratio 
Run for various Z' models and with a Z' mass of 1.5 TeV. The quoted statistical errors 
are those that can be expected for a luminosity of 100 fb^ 1 . 

4. Distinction between models and parameter determination 

Let us now discuss how well the different Z' models can be distinguished experimentally 
using the observables defined before: af[ x T, Ap B on- and off-peak, as well as R ua as 
obtained from the rapidity distribution. As a working hypothesis, a luminosity of 100 fb" 1 
and a Z' mass of 1.5 TeV will be assumed in the following. 

A precise knowledge of the cross section times the total width allows a first good 
distinction to be made between some models, as shown in the upper two plots of Fig. 6. 
It is not obvious how accurately absolute cross sections can be measured and interpreted 
at the LHC However, following the procedure outlined in [17], comparable reactions, in 
this case Z' and Z boson production, should be counted with respect to each other. The 
use of such ratio measurements should allow us to minimize systematic uncertainties, and 
an accuracy of ±1% might be achievable. As can be seen from the other plots in Fig. 6, 
the additional variables show a different sensitivity for the different couplings. 

For example, very similar cross sections are expected for the E 6 Z' models with cos (3 ~ 
±1 and for left-right models with a LR < 1.3. However, these two different models show a 
very different behaviour for on- and especially off-peak asymmetries and for the couplings 
to up-type and down-type quarks. Obviously, the maximum sensitivity can be obtained 
by using all observables together. Having said this, one also needs to point out that some 
ambiguities between the different models remain, even after a complete analysis of 100 
fb _1 of data. 

Assuming that a particular model has been selected, one would like to know how well 
the parameter (s), such as cos j3 or am, can be constrained. In the case of the E 6 model for 
instance, one finds that cos (5 cannot always be determined unambiguously. Very similar 
results can be expected for different observables but using very different values for cos (5. 
Again, the combination of the various measurements helps to reduce some ambiguities. 
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Figure 6: Variation of aff ■ T, A°p^ peak , A°p^ pcak and the ratio R q q as a function of the E 6 
model parameter cos/5 (left) and the LR-model parameter Qlr- The points corresponding 
to the particular Z' models are also shown. 
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If the Z' mass is increased, the number of events decreases drastically and the differ- 
ences between the models start to become covered within the statistical fluctuations. For 
the assumed luminosity of 100 fb -1 , we could still distinguish a Z' x from a Z' LR over a 
large parameter range; the Ap B measurements provide some statistical significance up to 
Mz> = 2-2.5 TeV. On the contrary, a Z' v could be differentiated from a Z 1 ^ only up to a 
Z' mass of at most 2 TeV as, in that case, the dependence of Ap B is almost identical in 
the two models. 

In summary, a realistic simulation of the study of the properties of Z' bosons origi- 
nating from various theoretical models has been performed for the LHC We have shown 
that, in addition to the Z' production cross section times total decay width, the measure- 
ment of the forward-backward lepton charge asymmetry, both on the Z' peak and in the 
interference region, provide complementary information. We have also shown that a fit 
of the rapidity distribution can provide a sensitivity to the Z' couplings to up-type and 
down-type quarks. The combination of all these observables would allow us to discrim- 
inate between Z' bosons of different models or classes of models for masses up to 2-2.5 
TeV, if a luminosity of 100 fb -1 is collected. 

Acknowledgements: We thank J. Hewett, G. Polesello and T. Rizzo for helpfull dis- 
cussions during the Les Houches Workshop. 

References 

[1] The LEP Electroweak Working Group and the SLD Heavy Flavour Group, Note 
LEPEWWG/2003-01; http : //lepewwg . web . cern . ch/LEPEWWG. 

[2] See for instance, M. Cvetic et al., Phys. Rev. D56 (1997) 2861. 

[3] J. Hewett and T. Rizzo, Phys. Rep. 183 (1989) 193; A. Leike, Phys. Rep. 317 (1999) 
143; M. Cvetic and P. Langacker, hep-ph/9707451. 

[4] See for instance, C. Hill and E. Simmons, hep-ph/0203079. 

[5] N. Arkani-Hamed et al., JHEP 0208 (2002) 021. For the phenomenological aspects, 
see: T. Han et al., Phys. Rev. D67 (2003) 095004; J. Hewett, F. Petriello and T. 
Rizzo, hep-ph/0211218; C. Csaki et al., hep-ph/0303236. 

[6] R. Casalbuoni et al., Phys. Lett. B460 (1999) 135; J. Rosner, Phys. Rev. D61 (2000) 
016006; J. Erler and P. Langacker, Phys. Rev. Lett. 84 (2000) 212. 

[7] For an account, see J. Erler and P. Langacker, Phys. Lett. B456 (1999) 69; see also 
for an updated analysis of LEP1 limits, F. Richard, hep-ph/0303107. 



13 



[8] M. Cvetic, S. Godfrey et al., DPF study in "Electroweak Symmetry Breaking and 
Beyond the Standard Model", hep-ph/9504216; T.G. Rizzo, hep-ph/9612440; S. 
Godfrey, in APS/DPF/DPB Snowmass Study, hep-ph/0201092 and hep-ph/0201093. 

[9] R. Robinett and J. Rosner, Phys. Rev. D25 (1982) 3036; (E) ibid. D27 (1983) 679. 

[10] Proc. of the ECFA Large Hadron Collider Worshop, Aachen (Germany) 1990, Reports 
CERN 90-10 and ECFA 90-133. 

[11] CMS Collaboration, G.L. Bayatian et al, Technical Proposal, CERN/LHCC 94-38, 
LHCC/P1, 15 Dec. 1994. 

[12] ATLAS, Letter of Intent, CERN/LHCC 92-3 and ATLAS Collaboration, W. Arm- 
strong et al.; Technical Proposal, CERN/LHCC 94-43, LHCC/P2, Dec. 1994. 

[13] A. Djouadi et al., Z. Phys. C56 (1992) 289. 

[14] P. Langacker, R. Robinett and J. Rosner, Phys. Rev. D30 (1984) 1470; V. Barger et 
al., Phys. Rev. D35 (1987) 2893; J. Rosner, Phys. Rev. D35 (1987) 2244. 

[15] M. Dittmar, Phys. Rev. D55 (1997) 161. 

[16] F. del Aguila, M. Cvetic and P. Langacker, Phys. Rev. D48 (1993) R969. 

[17] M. Dittmar, F. Pauss and D. Zurcher, Phys. Rev. D56 (1997) 7284. 

[18] T. Sjostrand et al, Comput. Phys. Commun. 135 (2001) 238. 

[19] G. Azuelos and G. Polesello, hep-ph/0204031; see also, T. Rizzo, hep-ph/0305077. 

[20] J.Hcwett and T. Rizzo in Ref. [3]. See also, A. Djouadi, Z.Phys.C63 (1994) 317; 
A. Djouadi and G. Azuelos, Z.Phys.C63 (1994) 327. 

[21] H. L. Lai et al. [CTEQ Collaboration], Eur. Phys. J. C 12 (2000) 375. 

[22] G. Altarelli, R.Ellis and G. Martinelli, Nucl. Phys. B 143 (1978) 521 and B146 (1978) 
544; J. Kubar-Andre and F. Paige, Phys. Rev. D19 (1979) 221. 



14 



